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J)  ABSTRACT 

Herein,  we  propose  a  mathematical  theory  of  thermo-viscoplasticity 
which  is  a  synthesis  of  experimentally  observed  material  behavior  on  one 
hand,  and  the  concepts  of  irreversible  thermodynamics  on  the  other. 

The  underlying  principle  is  that  the  history  of  deformation  is  de¬ 
fined  in  terms  of  a  "time  scale"  which  is  not  measured  by  a  clock,  but 
is  in  itself  a  property  of  the  material  at  hand. 

The  theory  is  unifying  in  the  sense  that  theories  of  plasticity, 

•  viscoelasticity  and  elasticity  can  be  obtained  from  it  as  special  cases 
by  imposing  suitable  constraints  on  the  material  parameters  involved, 
furthermore,  it  does  not  roaJce  use  of  the  idea  of  a  yield  surface. 
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Herein,  we  propose  a  mathematical  theory  of  therno-viscoplasticity 
which  is  a  synthesis  of  experimentally  observed  material  behavior  on  one 
hand,  and  the  concepts  of  irreversible  thermodynamics  on  the  other. 

The  underlying  principle  is  that  the  history  of  deformation  is  de¬ 
fined  in  terms  of  a  "time  scale"  which  is  not  sieasured  by  a  clock,  but 
is  in  itself  a  property  of  the  material  at  hand. 

The  theory  is  unifying  in  the  sense  that  theories  of  plasticity, 
viscoelasticity  and  elasticity  can  be  obtained  from  it  as  special  cases 
by  imposing  suitable  constraints  on  the  material  parameters  Involved j 
furthermore,  it  does  not  make  use  of  the  idea  of  a  yield  surface. 


Section  1.  Lndochronic  Theory  of  Viscoplasticity 

In  current  theories  of  plasticity,  to  explain  the  observed  discontin¬ 
uities  in  material  behavior  upon  loading  beyond  the  "yield  point"  and  upon 
unloading,  one  has  to  introduce  the  concept  of  a  yield  surface  in  stress 
space  as  well  as  a  "loading  function"  to  distinguish  between  loading  and  un¬ 
loading.  Similarly,  in  the  case  of  viscoplasticity ,  the  existence  of  a 
static  stress-strain  relation  and  a  yield  surface  are  assumed  and  the  stress 
increment,  with  respect  to  the  static  value,  is  related  to  the  strain  rate, 
or  more  generally  to  the  strain  history,  by  a  constitutive  equation. 

However,  the  fact  that  the  phenomenon  of  yield  is  usually  a  gradual 
transition  from  a  linear  to  a  non-linear  stress-strain  response,  makes  it 
difficult  to  say  precisely  where  yield  has  occurred,  to  the  extent  that 
different  definitions  of  yield  are  used  for  this  purpose.  Three  such  Defin¬ 
itions,  for  instance,  are  (a)  the  deviation  from  linearity  in  the  relation 
between  some  measure  of  strain  and  stress  (b)  the  interaction  of  the 
initial  part  of  a  stress  strain  curve  and  the  backward  linear  extrapolation 
of  the  "plastic"  part  of  the  curve  and  (c)  a  value  of  "proof"  stress 
corresponding  to  an  arbitrarily  defined  value  of  "oroof"  strain. 

Though,  from  an  engi?ieering  viewpoint,  the  initial  yield  surface  is  not 
overly  influenced  by  the  definition  of  yield,  it  has  been  found  experimentally 
that  subsequent  yield  surfaces  of  a  strain  hardening  material  are  influenced 
by  the  definition  of  yield  to  an  extraordinary  degree.  (See  Appendix  I).  If 
we  insist  that  the  increment  in  plastic  :s crain  is  to  be  normal  to  the  yield 
surface,  then,  for  complex  stress  histories,  each  such  definition  will  give 
rj.se  to  a  different  plastic  strain  history.  Only  one  of  these  can  be  the 


-correct  one . 
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The  conceptual  difficulties  that  are  encountered  by  the  introduction  of 
the  yield  surface  are  completely  circumvented  by  our  theory  cf  plasticity 
which  is  developed  on  the  basis  of  the  observation  that  the  state  of  stress 
in  the  neighborhood  of  a  point  in  a  plastic  material  depends  on  the  set  of 
all  previous  states  of  deformation  of  that  neighborhood,  but  it  does  not 
depend  on  the  rapidity  at  which  such  deformation  states  have  succeeded  one 
another* . 

The  independence  of  stress  of  the  rapidity  of  succession  of  deformation 
states  is  achieved  by  introducing  a  time  scale  f,  which  is  independent  of  t, 
the  external  time  measured  by  a  clock,  but  which  is  intrinsically  related  to 
the  deformation  f  the  material. 

Of  course  there  are  many  ways  of  introducing  such  a  time  scale.  However, 
it  appears  almost  mandatory  that  4  should  be  a  monotonicallv  increasing  func¬ 
tion  of  deformation,  otherwise  two  different  states  of  deformation  could  exist 
"simultaneously"  i.e.  for  the  same  value  of  f, .  Furthermore,  a  positive  rate 
of  change  ,  of  the  internal  energy  density  c  with  respect  to  4  could  not  be 
interpreter  unambiguously  as  a  process  of  increasing  c,  if  d4  could  be  negative. 

A  logical  definition**  for  4  is  then  given  by  the  relation 


*In  the  present  Section  and  in  subsequent  Sections  (with  the  exception  of 
Section  2 )  we  shall  assume  that  mechanical  changes  take  place  in  a  constant 
temperature  environment,  such  as  an  isothermal  atmosphere.  The  thermal 
changes  in  the  material  will,  therefore,  be  mechanically  induced  and,  in 
general,  will  remain  small.  Conversely,  only  thermal  changes  of  this  nature 
will  be  considered  in  this  paper. 

^Alternative  but  less  general  definitions  have  appeared  in  the  litera¬ 
ture.  For  instance,  Ilyushin^ 1^  and  later  Rivlin^2^  defined  a  "time"  s  by 

the  relation  ds2  =  dC..dC...  However,  we  have  found  that  this  definition  is 

n  -3 

too  narrow  to  describe,  quantitatively,  material  behavior  in  the  plastic 
range  as  wilx  be  discussed  later.  The  effect  of  temperature  on  t,  will  at 
this  time,  be  considered  sufficiently  small  to  be  negligible.  For  a  more 
vague  allusion  to  this  possibility  see  auso,  Schapery  (3). 


3 


d‘2  PiiWdCiidCk£ 


(3.1) 


where  C..  is  the  Right  Cauchv-Green  tensor  and  P...  »  is  a  fourth  order  tensor 
rj  -  -  like 

which  could  depend  on  .  The  positive  definite  nature  of  d£2  requires 
that  be  positive  definite.  In  case  of  small  deformation 

d£2  =  p .  ,idc .  .de,  d 
i]  Vi 

where  is  the  small  deformation  strain  tensor  and  couiu  depend 


on 


e  .  .  . 

Actual  materials,  on  the  other  hand,  do,  in  general,  depend  on  the  his¬ 
tory  of  deformation  as  well  as  on  the  rapidity,  or  rate,  at  which  deformation 
states  succeed  one  another.  To  describe  materials  of  this  type  one  may  con¬ 
struct  a  theory  of  viscoplasticity  by  introducing  a  time  scale  1 ;  which  is  re¬ 
lated  to  the  external  time  t. 

It  appears  logical  to  define  ;  by  the  relationship 

d;2  =  a2dC2  +  Bzdt2 

where  a  and  B  are  scalar  material  parameters.  Henceforth  d C  will  he  called 

d  2 

an  "intrinsic  time  measure",  and  z(t) ,  such  that  >  0  (0<_c«»),  will  be 
called  an  "intrinsic  time  scale". 

In  our  theory,  the  stress  (among  cf^er  properties)  is  necessarily,  a 
functional  of  the  strain  history,  defined  with  respect  to  the  intrinsic  time 
scale,  the  latter  being  a  property  of  the  material  at  hand.  As  a  result  we 
have  called  our  theory  an  endochronic  theory  of  viscoplasti city . 

The  theory  will  now  be  developed  in  a  general  thermodynamic  framework  in 
Section  3.  Before  this  is  done,  however,  the  thermodynamic  foundations  are 
laid  in  Section  2. 


m»w  Tg  aaaaa 


aa — - 
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Section  2.  Thermodynamic  Foundations 


The  following  are  the  fundamental  laws  of  thermodynamics,  which  apply 
to  all  continuous  media  irrespective  of  their  constitutive  properties.  (for 
materials  that  are  solid-like,  in  the  sense  that  they  have  a  memory  of 
their  initial  configuration,  it  is  more  convenient  to  express  these  laws  in 
the  material  coordinate  system  xx).  In  differential  form,  these  are  the 
first  law  of  thermodynamics , 


e  =  (p  /2p)  i  u..  -  h  .  t  C, 

O  1]  ,1 

the  rate  of  dissipation  inequality, 

6y  =  (p  /2p )  t11  C..  -  i  -  nO  >  0 
o  — 

and  the  heat  conduction  inequality 


(2.1) 


(2.2) 


-h10,.  >  0  .  (2.3) 

1  ~ 

The  symbols  in  eq.'s  (2.1),  (2.2)  and  (2.3)  have  the  following  meaning:  pQe 
is  the  internal  energy  per  unit  mass;  oq  and  p  are  the  initial  and  current 
mass  densities  respectively;  t1-1  is  the  stress  tensor  in  the  material  coor¬ 
dinate  system  x1  ;  is  the  right  Cauchy-Green  tensor,  h1  is  the  heat  flux 
vector  per  unit  undeformed  area  in  the  material  system;  p^Q  is  the  heat  supply 
per  unit  mass;  9  is  the  temperature,  y  the  irreversible  entropy  and  ^  and  n 
are  the  free  energy  and  entropy,  respectively,  per  unit  undefermed  volume 
finally  a  subscript  following  a  comma  denotes  differentiation  with  respect 
to  the  corresponding  material  coordinate.  A  dot  over  a  quantity  denotes 
material  derivative  with  respect  to  time.  To  avoid  repetitious  statements, 
henceforth  we  shall  refer  to  C. ^  as  the  "deformation". 


In  the  case  of  dissipative  materials  the  stress,  the  internal  energy 
and  entropy  densities  (and,  therefore,  the  free  energy  density)  of  a  material 
neighborhood  depend  on  the  entire  history  of  deformation  and  temperature  of 
tnat  neighborhood. 

In  the  theory  o~  irreversible  thermodynamics  the  effects  of  history  are 

taken  into  account  by  specifying  that  the  stress  and  free  energy  density  are 

functions  of  the  current  values  of  C . .  and  0  as  well  as  n  additional 

il 

independent  variables  q  ,  not  necessarily  observable,  called  "internal  vari- 

a 

ablcs",  These  mav  be  scalars  or  components  of  vectors  or  tensers  ir.  the 
materitil  frame;  whatever  their  geometric  nature  they  must  remain  invariant 
with  translation  and  rotation  of  the  spatial  system  to  satisfy  the  principle 
of  material  indifference.  Thus: 


*  =  *  (C^  ,  0  ,  qa) 


=  T13  (Ck£  ,  e  ,  qo> 

(4) 

It  has  been  shown  elsewhere  that 


(2.4) 

(2.5) 


ij  _  2o  3 i|i 

T  o  3C.  . 
o 


n  = 


3i|> 

36 


0Y  =  -  q  >0 
,qa  Ma  - 


(2.6) 

(2.7) 

(2.B) 


Furthermore  the  heat  flux  vector  h1  is  a  function  of  8,.  ,  6  ,  C^..  and 


q  i.e. 
a 


h  =  h  (6^  -  6  .  ck£  ,  qo) 


(2.9) 


subject  to  the  conditions: 


dered  in  the  next  Section. 
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Section  3.  Constitutive  Liquations  in  Viscoplasf'  city 


from  the  right  hand  side  of  eq.  (7.3)  and  the  fact  that  jj* 
dz/dc  >  0,  it  follows  t fiat 


dt 

>  0 


^>q  d: 
o 


>  0  ('i  not  summed) 


(3.1) 


a  .  ■ 

where  the  inequality  is  valid  unless  — —  =  0.  It  also  follows  from  inequality 

dq  "  ^  ^  dqa 

(3.1)  that  -r~—  ,  u  ,  C.  •  and  o  must  be  related  otherwise  and  coul^ 

dz  u  in  3q  dz 

J  u 

be  prescribed  independently  and  in  such  a  fashion,  that  the  inequality  would 
be  violated.  In  this  event  there  must  exist  a  set  of  relations 


ar  -  f«(cii  •  •  e) 


(3.2) 


for  all  a,  where  the  functions  f  are  material  functions. 

It  must  be  noted  that,  as  a  result  of  eo .  (3.2)  q^  are  indeed  functionals 
of  the  histories  of  deformation  and  temperature  with  respect,  however,  to  the 
intrinsic  time  suic  r.  which  is,  itself,  a  material  property. 

Thus,  at  least  formally,  the  constitutive  equations  of  the  endocnronic 
theory  of  viscoplasticity  are  now  complete  in  the  sense  that  given  the  mate¬ 
rial  functions  f  and  h1  then  for  some  specified  deformation  and  temperature 
a 

histories,  q  are  found  from  eq .  (3,2)  and  thus  t1"1  and  n  are  found  from  eq's. 
(2.6)  and  (2.7)  respectively;  similarly  h1  the  heat  flux  vector  is  determined 
from  eq .  (2.9). 

Ideally,  one  would  like  to  know  the  thermomechanical  three-dimensional 
response  of  a  material  over  the  whole  spectrum  of  mechanical  and  thermal  con¬ 
ditions,  i.e.,  under  all  variations  in  strain,  strain  rate  (or  more  generally, 
history  of  strain)  and  temperature.  However,  such  a  task  would  be  a  momentous  * 
if  not  an  impossible,  undertaking;  the  experimental  evaluation  of  the  material 
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functions  involved,  under  wide  conditions  of  strain  and  temperature,  would  be 
impract i cal . 

Fortunately  the  domain  of  specification  of  design  conditions  is  usually 
limited  in  some  way;  for  instance  usually,  (a)  large  changes  of  temperature, 
fast  rates  of  loading,  but  small  strains  are  prescribed;  or  (b)  small  changes 
in  temperature  and  small  rates  of  loading  but  large  strains  and/or  displace¬ 
ments  prevail.  More  extreme  mechanical  as  well  as  thermal  conditions  are 


It  is  reasonable  to  expect  that  material  behavior  would  be  easier  to 
describe  mathematically  over  a  narrower  domain  of  environmental  conditions, 
where  the  applicability  or  "correctness"  of  such  mathematical  formulation 
would  be  easier  to  check  experimentally. 

In  what  follows  we  shall  consider  situations  in  which  the  strain  in  a 
material  region  R  as  well  as  the  temperature  changes  relative  to  a  uniform 
reference  temperature  are  "small".  To  make  the  above  statement  more  pre¬ 
cise  let  c.  .(z')  denote  the  history  of  the  strain  tensor*,  for  z  <  z'  <  z 
13  o  -  — 

where  z  is  some  initial  intrinsic  time.  Set 
o 


|e..(z')|  |  Mc.^(z')  €ij(z^  ' 


(3.3) 


and  let  the  supremum  of  | | z') | | be  A. 

Similarly  let  SKz')  be  the  history  of  the  temperature  increment  relative 
to  the  reference  temperature  and  let  the  supremum  of  |9(z')|  be  6.  The 
notion  of  smallness  is  made  precise  by  stipulating  that  A<<1  and  6<<1. 

Thus,  formally 


e  +  S 
0 


(3.4) 


(C.  .  -  6.  .) 

-  j  y  j 


9 


n  =  no  +x 


(3.5) 


e.  ,(z')| j  =  A  ,  |0(z')|  =  « 

1  i]  ' sup  sup 


(3.6) 


To  complete  the  formalism  let  x  be  the  entropy  change  relative  and  a 


reference  uniform  entropy  nQ  ,  and  let  o„  denote  the  stress  tensor.  The  re¬ 
ference  state  is  defined  by  the  condition  that  =  0,  =  0  =  x  =  0,  qQ  =  0. 

Under  these  conditions,  eq.’s  (2.4)  through  (2.9)  and  eq.  (2.10b)  become. 


* =  +  (eij  *  J  >  (’a) 

(3.7) 

0.,  =  |t 

ID  3e .  . 

ID 

(3.8) 

X 

II 

1 

(3.9) 

er  =  -  1^  q  >0 
o  3q  — 

(3.10) 

h.  =  k..  0  ,. 
l  1)  D 

(3.11) 

h.  ,.  <  0 

(3.12) 

where  is  the  thermal 
Finally  eq.  (2.11), 


conductivity  tensor, 
becomes 


+  Q 


(3.13) 


It  is  shown  in  Appendix  II  that  |aol  and  Iq^J  may  stay  small  in  the  sense 
that  given  two  positive  numbers  6^  and  6°  ,  however  small,  A  and  6  can  be  cho¬ 
sen  small  enough  such  that  |qal  <  6^  and  |qQ|  <  6° 

At  this  stage  one  may  obtain  the  corresponding  equation  for  qfl  by 
linearizing  eq.  (3.2).  However  it  is  physically  more  meaningful  and,  a  pos¬ 
teriori,  more  rewarding  to  examine  more  closely  the  rate  of  change  of 
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irreversible  entropy  y 
From  eq  .  (3.10), 


e  p.  =  _  ^  >  o 

o  dz  3q  dz  — 

a 


(3.14) 


It  follows  from  eq.'(s)  (3.14)  and  (3.?)  that  ^  may  be  expressed  as  a  func- 

dq  j 

tion  of  -t  -a  ,  c.  .  and  ,  subiect  to  the  condition  that  t —  =  0  whenever 
,  dz  il  dz 

dq  d 

-  0  for  all  a  .  Thus,  if  we  expand  G  in  a  Tavlor  series  and  ignore 

dz  *  o  dz 

terms  of  order  higher  than  0(62)*  and  observe  the  inequality  (3.10),  thereby 

a 

eliminating  the  linear  terms  in  the  expansion,  then 


,  dq  do. 

e  dy.  -  b  8 

o  dz  afi  dz  dz 


(3.15) 


Eq.'s  (3.14)  and  (3.15)  are  simultaneously  satisfied  if 


♦  b  „  o 

3q  a6  dz 

a 


(3.16) 


With  eq.  (3.16)  the  constitutive  description  of  a  viscoplastic  material  is 
now  complete. 

Explicit  Constitutive  Equations 

Explicit  constitutive  equations  for  viscoplastic  materials  under  condi¬ 
tions  of  small  strain  and  small  changes  in  temperature  are  obtained  bv  ex¬ 
panding  4>  in  eq.  (3.7)  in  Taylor  Series  and  omitting  terms  of  order  higher 
than  0(62);  *  linear  terms  must  vanish  to  satisfy  the  initial  conditions. 

Before  the  expansion  is  carried  out,  however,  it  appears  desirable  to 
regard  q  not  as  scalars  but  components  of  second  order  tensors.  This,  as 

Q 

will  be  shown,  obviates  certain  difficulties  which  arise  with  the  represen¬ 
tation  of  fourth  order  tensors.  For  instance,  in  Fef .  (**),  we  were  faced 


*6  is  the  largest  of  6^  and  6^.  Also,  5  is  the  largest  of  5,  A  and  6.. 
a  1  i 
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with  having  to  assume,  without  proof. 


that  a  fourth  order  tensor 


Ci  jkf 


such 


that , 


*"ijk£  ('jik£  ^ij£.k  ^kiij 


i3  given  by  the  series 


C.  ,  =  E 
1  jk-t  a 


H  .  .  n 

k4q 


0.17) 


where  a.,  are  second  order  symmetric  tensors  and  a  are  scalars.  Problems 
l  ja  a 

such  as  this  are  obviated  by  giving  the  internal  variables  a  tensorial  char¬ 
acter,  Thus  the  free  energy  density  and  other  thermodynamic  quantities  are 
now  functions  of  c . .  ,3  and  n  internal  variables  q? .  (a  =  1 ,  2  .  .  .  .  n) , 

l]  — 

where  q?..  are  symmetric  second  order  tensors,  with  respect  to  the  material 
system  .  In  this  notation,  eq.’s  (3.13)  through  (3.16)  now  read. 


h.  .  =  e 

i,i  o  nS 


3£ 

3q 


•a  i 
q.  .  +  Q 

a 

ij 


e  ^ii>  o 

o  dz  .  a  dz  — 

3qij 


(3.18) 


(3.19) 


.  dy  ...  a  ,  a  .a 

0  dz  a  ijk£  -r^ii  -r^k£ 

J  dz  '  dz 


(3.20)* 


and 


h«  dqk^  _ 
.  a  ijk£  dz 
3qij 


(a  not  summed) 


(3.21) 
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*Expansions  of  the  type  b... 
the  above  form.  See  Ref.  5. 
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.  0  -z - -r -  and  A .  „  q.  .  q.  »  reduce  to 

i3ki  dz  dz  i]kc  ^k£ 
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Furthermore,  in  view  of  my  previous  discussion. 


*  =  ^ijk*  Cij  ek£  *  Bijk£  cij  % i  +£Cijk*  qij  V 


+  D.  .  $  e.  .  +  E?.Sq®.  +  *sF  5  2 
iD  J-D  i-D  ID 


(3.22) 


Though,  in  principle,  eq.'s  (3.8),  (3.9),  (3.18),  (3.21)  and  (3.22)  are 
sufficient  for  the  derivation  of  explicit  constitutive  equations,  we  shall 
obtain  these  only  for  isotropic  material,  so  as  to  keep  the  algebra  at  a 
minimum.  For  such  materials 


A...,  *  A,  6..  6,.  +  A„  6..  6./, 
ljk-t  1  id  kt  2  lk 


B®. ,  .  =  B®  6 .  .  6.  .  +  B®  6..  6.„ 
lDk 1  1  i]  k£  2  lk  d Z 
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D.  .  =  D  <5.  . 

ID  ID 


E?  .  »  E°  6.  . 

ID  ID 


k.  .  =  k  6. . 

ID  ID 


b®  „  =  b®  <5.  .  +  b® 

ijkt  1  i]  kc  2  ik  34 


(3.23  a-f) 


It  is  worth  noting  that  here  we  consider  materials  which  are  "stable" 

in  the  sense  that  straining  of  the  reference  configuration  under  isothermal 

conditions  will  increase  the  free  energy  density  ^  .  Thus  A^^  and 

are  positive  definite.  As  a  consequence  A^  ,  A^  ,  C®  and  C®  are  all  positive. 

Omitting  superfluous  algebra,  the  coupled  thermomechanical  constitutive 

°kk 

equations  take  the  folVWng  form  in  terms  of  the  hydrostatic  stress  o  *  , 

the  deviatoric  stress  te  -cor  ck..  ,  the  increment  in  temperature  3  ,  the 


hydrostatic  strain  c,  ,  ,  the  deviatoric  strain  tensor  e..  and  the 

1  kk  il 

entropy  increment  x  >  in  terms  of  the  above  notation: 


s  .  .  =  2 
i] 


p  ( z-z'  )  j  7-*-  dz' 


(3.24) 


"  32kk  lz  as 

K(z-z')  .  ■  dz*  +  D(z-z')  dz' 

1  z  2  ^  z 

0  0 


(3.25) 


-X  5  =  [  D(z-z')  kk  dz'  +  F(z-z')  — ,  dz'  <3.26) 

j  U  _ 


where , 


B2  B2  B2  B2  "paz 

2  u(z>  =  (A2  -  l  )  H(z)  ♦  l  -H-J.  e  a 


a  c; 


a  c: 


(3.27) 


Ba  Ba  Ba  Ba  -X  z 

K  (z)  =  (A  -  7  -2-2-  )  H(z)  ♦  i  — — -  e  a 

1  °  L  Ca  a  Ca 

o  o 

BE  B  E  -X  z 

D(2)  =  (d  -  T  -2—1  )  H(  2 )  ♦  l  -2_“  e  ° 

a  C°  a  C° 

o  o 

ra  ra  a  a  -X  z 

r(z)  =  (F  -  l  )  H(  z)  e  a 

a  Co  a  Cc 


(3.28) 


(3.29) 


(3.30) 


A  s  1/3  ( 3A.  +  A„)  ,  B  =  1/3  (3B,  ♦  B„)  ,  C  =  1/3  (3C.  ♦  C. )  (3.31) 

0  I/O  1/0  1/ 


.  S  x  _  fo 
a  ”  .a  ’  a  ~  .  a 
b2  bo 


(3.32) 


»:  ■  3b°  *  b2 


(3.33) 


The  heat  conduction  equation  is  similarly  found  to  be: 


dt  .  q  _  3  fZ  r  .  , .  39  ,  ,  .  3  ( 
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_  ,  r,  a  a  a  r,  a 

Q  +  Ib1  q..  q..  ♦  Ib2  qtj  q£j 


where 

\ 


C'f-)  =  -  e  f(z) 
v  o 


(3.35) 


and  *  over  a  quantity  implies  differentiation  with  respect  to  z.  The 

lo«er  limit  z  denotes  the  intrinsic  time  of  the  reference  state. 

i  ® 
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Section  4.  Endochronic  Theory  of  Plasticity  and  its  Relation  to  Present  Theories 


Our  theory  of  plasticity,  which  is  a  rate  independent  endochronic  theory, 
is  obtained  by  replacing  the  time  measure  dt  by  d£.  The  time  scale  now  becomes 
z(€),  but  the  form  of  the  constitutive  equations  remains  unaltered.  In  parti¬ 
cular,  the  "linear"  form  of  our  theory  is  obtained  by  setting 


2 


=  Pijkt  dEij  dck£ 


(4.1) 


where  Pijj^  is  a  positive  definite  fourth  order  material  tensor.  We  repeat 
the  constitutive  equations  of  the  linear  theory,  in  the  particular  case  when 
the  deformation  is  isothermal  so  that  a  comparison  may  be  made  with  current 
theories.  When  9=0,  then  eq.'s  (3.24)  and  (3.25)  become. 


(4.2) 

(4.3) 


where  z  ■  z(t). 

If  the  material  behaves  elastically  under  pressure  (so-called  plasti¬ 
cally  incompressible)  then  K(z )  is  a  constant  and  in  this  case 


°kk  =  3Kckk 


(4.4) 


Whereas , 


s .  .  =  2 
13 


u  (z-z*  )  de^iz' ) 


(4.5) 


Let  now  u  (z )  consist  of  a  single  exponential  term  i.e. 


4  (z  )  =  u  e 


(4.6) 
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In  this  event 


3ij  ■  J“o 


-a(z-z' ) 


de^  (s') 


(4.7) 


The  integral  eq.  (4.7)  is  reducible  to  the  differential  equation 


de;  j  =  dz  s.  .  ♦  — rr  ds. 


'  27o  az  sij  T  TTo  as“ij 


(4.8) 


Now,  ds„  may  be  identified  as  the  "elastic"  component  of  deviatoric 
strain  of  classical  plasticity.  If  one  follows  the  traditional  definition 
of  "plastic  strain"  de^P  given  below,  i.e.. 


d«fjP  =  de  •  j  -  de®.. 


then  in  view  of  eq.  (48) 


deijP  *2irodzaii 


(4.9) 


(4.10) 


But  these  are  the  Prandtl-Reuss  relations.  Hence  our  present  theory  contains 
these  relations  as  a  special  case.  Where  then  does  it  differ  from  this  theory? 
It  does  in  the  interpretation  of  the  proportionality  coefficient  dz.  In  the 
Prandtl-Reuss  theory  may  be  positive  negative  or  zero,  and  in  fact,  z  has 
been  identified  with  the  yield  surface,  i.e.,  plastic  action  is  assumed  to 
occur  when  dz  %  o ,  where 


z  =  Z(s..) 


but  that  de^j*3  is  zero  whenever 


dz  <  o 


(4.11) 


(4.12) 


In  the  present  theory  dz  is  always  positive  if  the  material  is  deforming  (it 
is  zero  only  when  deformation  does  not  take  place).  Thus  always , 
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dz  >  O 


(4.13) 


furthermore  dz  is  not  given  by  eq.  (4.11)  i.e.  it  is  rot  related  to  somo 
yield  surface  but  its  definition  is  entirely  kinematic.  Thus,  no  yield 
phenomenon  or  surface  are  postulated  here.  One  obtains  the  stress  response 
by  merely  monitoring  the  history  of  strain. 

Also  the  theory  admits  a  further  generality  since  u(z)  need  not  consist 
of  a  single  exponential  term. 

For  instance  u  (z)  may  be  of  the  form 

u  ( z  )  =u  +u  ^e~aZ  .  (4.14) 

In  this  case,  however,  the  differential  form  of  eq .  (4.7)  becomes: 


2(u.  ♦u,)de..+2uftae..ds  =  ds..+cs,,dz.  (4.15) 
Oil)  0  1^  i} 

The  shear  modulus  u,  at  z  =  0,  (initial  modulus),  is  (u  +u  ^ ) .  The 
"plastic"  components  of  the  deviatoric  shear  strain  tensor  are  given  from 
eq.  (4.9),  i.e., 


d*ii  * 


2  u 


d  z 
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-  2  p 


0  Bij’ 


(4.16) 


Note  that  eq .  (4,16)  does  not  satisfy  the  Frandtl-Reuss  relations,  which  are 
also  violated  if  one  adds  more  exponential  terms  to  the  right  hand  side  of 
eq.  (4.14).  In  fact  these  relations  win  be  satisfied  if  and  only  if  u  (z  ) 
is  given  by  eq.  (4.6),  i.e.,  u  is  represented  bv  a  single  exponential  term 
only.  This  situation  is  not  particularly  disturbing.  Peters  et  Als^**^ 
carrie.d  out  experiments  on  thin  walled  14S-T4  aluminum  alloy  cylinders  by 
loading  these  in  combined  compression  and  torsion  and  found  that  the  Frandtl- 
Feuss  relations  were  not  satisfied,  for  this  particular  metal. 


Conclus ions 


A  theory  has  been  presented  here,  the  scope  of  which  is  wide  enough  to 
allow  a  rational  phenomenological  description  of  mechanical  behavior  of 
materials  under  various  histories  of  strain  and  temperature.  In  particular, 
the  viscopiastic  behavior  of  materials  is  formulated  mathematically,  without 
recourse  to  the  dichotomy  of  the  deformation  history  in  plastic  and  elastic 
parts  and  without  the  necessity  of  introducing  discontinuities  in  material 
behavior,  such  as  yield  surfaces. 

The  theory  merely  asserts  that,  to  every  history  of  deformation  gradient 
and  temperature  of  a  neighborhood  there  corresponds  a  unique  state  of  stress 
in  that  neighborhood.  An  entirely  novel  feature  of  the  theofy  is  that  these 
histories  are  defined  with  respect  to  a  time  scale,  which  itself  is  a  mate¬ 
rial  property. 

In  this  paper,  we  have  merely  presented  the  framework  of  the  theory 
without  actually  evaluating  the  material  functions  involved,  through  the  use 
of  experimental  data.  This,  however,  will  be  done  in  Part  II  of  thi3  paper, 
where  it  will  be  shown  that  the  theory  describes  experimentally  observed 
plastic  behavior  of  metals  with  a  remarkable  degree  of  accuracy. 


Stipulation 


This  manuscript  is  submitted  for  publication  with  the  understanding  that 
the  United  States  Government  is  authorized  to  reproduce  and  distribute  re¬ 
prints  for  governmental  purposes. 


The  following  is  a  short  account  of  the  experimental  work  on  (a)  the 
effect  of  the  definition  of  yield  on  the  shape  of  the  yield  surface  and  (b) 
of  the  work  on  viscoplasticity.  The  references  given  are  by  no  means  exhaus¬ 
tive  and  the  author  wishes  to  apologize  to  people  of  whose  work  he  is  not 
currently  aware. 

In  Ref.'s  6  and  7,  aluminum  alloy  tubes**  were  subjected  to  shear  pre- 
strain  by  twisting  well  into  the  plastic  region  by  a  nredeterminod  amount. 

The  yield  surface  corresponding  to  this  degree  of  prestrain  was  established 
by  loading  the  tubes  in  combined  tension  and  torsion. 

In  Ref.  G,  Haghdi  found  that  subsequent  yield  surfaces  distorted  in  the 
direction  of  the  shear  axis  with  a  pronounced  Bauschinger  effect  in  shear 
but.  there  was  no  effect  on  the  yield  stress  in  tension  (i.e.  the  yield  locus 
did  not  change  in  the  vicinity  of  zero  shear  stress). 

On  the  other  hand,  in  Ref.  7,  Ivey  observed  that  the  yield  surface,  in 
addition  to  distortion,  underwent  a  large  amount  of  translation  ir.  the  direc¬ 
tion  of  the  shear  axis,  so  that  for  large  prestrains  the  origin  of  the  stress 
space  was  outside  the  yield  surface.  However,  ’  .  was  in  agreement  with 
h'aghdi  in  that  the  presence  of  shear  prestrain  did  not  affect  the  yield  stress 
in  tension.  Both  authors  used  deviation  from  linearity  in  a  stress  strain 

diagram  as  a  definition  of  vicld. 

(  8 ) 

Mair  and  Pugh  to  check  the  absence  of  "cross  effect",  carried  out 
their  own  experiments  on  copper  with  a  high  degree  of  isotropy.  However  they 
used  a  different  definition  of  yield,  this  being  the  point  of  intersection 
of  the  initial  straight  part  of  the  stress-strain  curve  with  a  backward  linear 
extrapolation  of  the  "plastic"  part  of  the  stress-strain  curve. 
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Their  results  varied  significantly  from  those  of  Ivey  and  Naghdi.  They 
found  that  expansion  and  distortion  of  the  initial  locus  took  place  with  a 
strong  cross  effect  between  shear  and  tension.  Also  a  pronounced  Bauschinger 
effect  in  torsion  was  found  with  large  initial  positive  pretorsion.  These 
authors  also  observed  pronounced  "plastic"  unloading  in  shear. 

(9) 

The  results  of  Szczepinski  and  Miastkowski  tend  to  confirni  the  find- 
( 8 ) 

ings  of  Mair  and  Pugh  .  Their  results,  moreover,  were  significant  in  other 
respects.  Specifically,  using  the  proof  strain  to  define  yield,  they  studied 
aluminum  alloy  sheets  under  biaxial  tension  with  the  intention  of  finding 
the  effect  of  prestrain  on  the  shape  of  the  yield  surfaces.  They  observed, 
migration,  distortion,  expansion  and  sometimes  rotation  of  the  initial  yield 
locus. 

Similar  conclusions*  can  be  drawn  from  Szczepinski 's  paper^10\  as  well 

( 17  ) 

as  Miastkowski  and  Szczepinski ' s  ,  in  which  tubular  brass  specimens  were 
subjected  to  combined  axi3l  and  circumferential  stress. 

Initial  and  subsequent  yield  loci  were  plotted  when  yield  was  defined 
(a)  as  departure  from  linearity  or  (b)  when  it  was  set  to  correspond  to  a 
certain  proof  strain.  In  particular,  when  definition  (a)  was  used,  subse¬ 
quent  yield  loci  did  not  contain  the  initial  locus,  but  when  (b)  was  used, 
with  proof  strain  set  at  0.5%  subsequent  loci  contained  the  initial  locus. 

Attempts  to  describe  tne  change  of  the  yield  locus  with  prestrain,  by 

-  (13) 

simple  models  have  not  proved  satisfactory .  Batdorf  and  Budianski  , 
suggested  that  after  prestrain,  the  yield  locus  is  the  minimum  surface 


i  ! 
I  i 

5  ; 


*  In  this  connection,  see  also  work  of  the  same  general  nature  by 
Bertsch  and  Findley^11'1  and  Ku  and  Bratt^1?\ 
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through  the  point  of  prestrain  and  the  initial  yield  locus.  This  model 

however  does  not  account  for  the  Bauschinger  effect.  The  kinematic  hardening 

( m ) 

rule**  proposed  by  Prager  ,  was  partially  successful,  m  so  far  as  it  can 
be  of  value  only  when  the  stress-strain  curve  of  a  material  in  simple  tension 
is  bilinear^^.  Otherwise  subsequent  shapes  of  yield  locus  must  be  defined 
in  terms  of  a  parameter  that  depends  on  the  .listory  of  strain^^  to  obtain 
realistic  unloading  behavior. 

A  more  realistic  model  is  the  one  by  Hodge^^  which  includes  transla¬ 
tion,  expansion  and  distorsion  of  the  yield  surface.  This  model  covers  all 
contingencies  but  does  not  include  the  history  of  stress  on  the  shape  and 
position  of  the  yield  locus. 

However,  every  definition  of  yield  gives  rise  to  a  different  yield  sur¬ 
face.  If  we  insist  that  the  increment  of  plastic  strain  is  tc  be  normal  to 
the  yield  surface,  then,  for  a  complex  but  specific  loading  history,  each 
such  definition  will  give  rise  to  a  different  plastic  strain  history.  Only 
one  of  these  can  be  the  correct  one. 

(15) 

So  it  appears  that  through  Eisenberg's  and  Phillip's  mathematical 
description  of  a  yield  surface  has  been  most  promising,  we  must  be  prepared 
to  question,  if  necessary,  whether  the  concept  of  yield  point  and  yield 
surface  are  the  only  way  by  which  plastic  effects  may  be  described,  especially 
in  view  of  the  fact  that  these  may  take  place  immediately  following  the  ini¬ 
tiation  of  deformation  of  material,  though  they  may  be  negligible  in  the 
region  of  small  strains.  This  would  agree  with  the  point  of  view  that  dis¬ 
locations  (and,  therefore,  plastic  behavior)  originate  immediately  upon 
initiation  of  the  loading. 


**  See  also  Ref.  (19). 


The  need  for  the  development  of  the  theory  of  viscoplasticity  arises 
from  the  recognition  of  the  strain  rate  sensitivity  of  metals  under  dynamic 
loading , 

The  difficulty  in  trying  to  synthesize  a  rational  "rate"  theory  from 
experimental  observations,  a  priori,  lies  in  the  fact  that  under  dynamic 
conditions  the  inertia  effects  are  significant.  In  the  absence  of  a  consti¬ 
tutive  theory,  these  effects  cannot  be  calculated;*  Therefore,  in  the  case 
of  dynamic  theories,  such  as  viscoplasticity,  theory  and  experiment  must 
advance  together. 

The  literature  abounds  with  data  on  the  subject  of  strain  rate  sensiti- 

(20-30)  (32) 

vity,  particularly  in  one  dimension.  Lindholn  carried  out  dynamic 

experiments  in  one  and  two  dimensions  in  an  attempt  to  generalize  results 
which  were  arrived  at,  by  consideration  of  thermally  activated  processes  and 
their  relation  to  dislocation  theory  in  metals.  See  Ref.'s  35-LO. 

An  early  attempt  at  a  theoretical  viscoplastic  constitutive  equation  in 

(31  32) 

one  dimension  is  due  to  Malvern.  ’  ~  This  equation  assumes  the  existence 
of  a  "static"  stress-strain  relation  and  then  relates  the  stress  increment, 
with  respect  the  static  value,  to  the  strain  rate. 

Modifications  and  generalizations  of  Malvern’s  equation  were  made  by 
Lubliner  who  included  a  limiting  maximum  stress-strain  curve,  and  by 
Perzyna^ 34  ^  ,  Perzyna  and  Wojno^  who  proposed  a  multiaxial  generalization 

*  Constant  strain  rate  experiments  would  appear  to  be  an  exception,  by 
being  less  susceptible  to  inertia  effects.  However,  Ref.  38  tends  to  negate 
this.  Long  specimens  give  different  responses  to  short  ones,  under  the  same 
conditions . 
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for  finite  strains,  assuming  the  additivity  of  the  elastic  and  plastic  strain 

( 39 ) 

components  and  by  Perzyna  who  used  concepts  of  internal  coordinates  and 
irreversible  thermodynamics  to  eliminate  the  above  assumptions  and  to  out  the 

theory  or.  sounder  foundations. 

.  ( 3° ) 

Though,  in  his  last  treatment,  Perzyna  abandoned  the  additivity  of 

plastic  and  elastic  strains,  he  still  retained  the  concepts  of  yield  stress 
(and  yield  surface)  and  the  hypothesis  of  a  datum  plastic  stress  strain  rela¬ 
tion,  with  respect  to  which  "strain  history"  is  to  be  related  to  the  "excess" 
stress  through  the  internal  coordinates. 

Our  theory  differs  from  Perzyna' s  theory  in  this  respect. 

We  close  by  mentioning  that  ,  with  the  exception  of  the  papers  by 
(  39 ) 

Perzyna  ,  only  a  moderate  research  effort  has  been  made  in  the  area  of 

coupling  hetween  a  viscoplastic  and  a  thermal  process.  However,  Chidister 

(25)  (22)  (37) 

and  Malvern  ,  Lindholm  and  Trozera,  Sherbv  and  Dorn  ,  considered 

the  effect  of  a  change  in  uniform  temperature  on  viscoplastic  behavior,  with 

a  view  to  confirming  some  results  of  the  dislocation  theory. 
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Appendix  II 

The  form  of  eq.  (3.22)  for  isotropic  materials,  according  to  eq.  *s  (3.22 
a-e),  is 

A2  ‘ii'ij  *  0-j 

’  '  «!i  <4*1  *  -  <$  ^  *■»<, 

+  h  F  9  2 


(A. 2.1) 


As  a  result  eq. 's  (3.8)  and  (3.9)  yield: 


°ij  *  A1  6ij  ekk  *  Vij  *  Bi  \k  ♦ 


+  D^5ij  sunune<1) 


X  =  |J  =  Deii  +  Eaq®.  +  FS  (a 


sunned) 


(A. 2. 2) 
(A. 2. 3) 


On  the  other  hand. 


^  *  B“l  Sk  \j  *  *2  *ij  *  C1  Sij  <& 
ij 

+  ^ 2  qij  +  9  (  not  summed) 


(A. 2. 3) 


Hence,  use  of  eq.  (3.21)  in  accordance  with  eq.  (3.23f)  yields  a  set  of  first 
order  differential  equations  for  qa^  ;  these  can  be  expressed  as  a  set  for 
qkk  and  •nother  for  the  deviatoric  part  of  q^ ,  which  we  denote  by  p“_  Thus 
in  the  notation  of  eq. 's  (3.31)  and  (3.33) 


a  dqkk 


8o  «kk  *  C”  >£k  ’  E“  *  b„ 


=  o 


(A. 2.4) 
(A. 2. 5) 
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In  both  eq.'s  (A. 2. 4)  and  (A. 2. 5)  a  is  not  summed.  It  follows  from  the 
abov(e  two  equations  that 


B®  tz 


a  tz 

4 ■  -  rl  | 

b  ;  % 


jJ.Jz-z') 


e.  ,  (z '  )dz  '  -  —  iX«(z"z,)-5(z')dz'  ( A .  2 .  ( 

kk  Ka  * 


b  *  z 
0  0 


p“  =  -  5pa^z_z  )ei.(z')dz'  (A. 

b2  zo 

where  \  and  p  are  given  by  eq.  (3.32). 
ol  a 

In  the  light  of  the  tensor ial  notation  that  we  have  adopted  for  the 


(A. 2. 7) 


internal  variables,  let 


,  M  | ! 

e .  .  s  4  ,  lit,, 

n  in 

sup  11  J 1 ' sup 


A1  •  iekki  =  Ao 


(A. 2. 8) 


where,  i|e„||  =  I c £ ^  cij'  *  etc-  Evidently  , 


(A. 2. 9) 


Then  as  a  result  of  eg.s  (A. 2. 7)  and  (A. 2. 8) 


<  1  2 1  A, 


(A. 2. 10) 


a  'Bcl  Ie“I  , 

qkk  1  “a  Ao  ♦  “S’  5 

0  0 


where  as  before  \9 


Also  from  eq «  ( A . 2 . 5 ) , 


(A. 2. 11) 


(bp2 


2  2.2 

=  (B  )  e.. 

a  11 


2  ,02, 2 

♦  (C  ) 
a 


(A. 2. 12) 


2  2 

t  2B  C  p. .  e. . 

a  a  'lj  i] 


However,  since 


P,  .  e. .  <  p. .  e. - 

ii  I  13 


(A. 2. 13) 
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it  follows  from  (A. 2. 12)  that 


Also  as  a  result  of  eq.  (A. 2. 4) 

,  a 

h°  |J^J  <2  A  +  2  E®  I  6 

o  |dz  |  —  oo  | 

At  this  point  we  order  our  internal  variables  as  shown, 

12  m  1  2  m 

pij  *  pij . pij  ;  qkk  ’  qkk  ’  ‘  '  qkk  * 


(A. 2. 14) 


(A. 2. 15) 


Let  be  a  typical  internal  variable.  Then,  whether  it  belongs  to  "he  p- 

group  or  the  q-group  above,  as  a  result  of  eq.'s  (A. 2. 10),  (A. 2. 11),  (A. 2. 14) 

and  (A. 2.14),  given  two  positive  members  4^  and  4  2 ,  however  small,  we  can 

choose  A  and  A,  (and  therefore  A)  and  4  such  that 
o  1 

dq 

qa  i  *1  and  dz^  -  62  * 
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